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I. INTRODUCTION 


Perhaps the most intuitive, and certainly one of the most useful 
characterizations of a random process is its power spectrum. Unfortunately, 
the power spectrum is also a somewhat elusive function in the sense that 
it cannot be determined unambiguously by observing a process over a finite 
time interval. So, the problem of how best to use available observations 
of the process to estimate its spectrum is one which has been of considerable 
interest to scientists and engineers for many years. 

Basic aspects of the spectral analysis problem will be briefly reviewed 
in this report, but the primary interest here will be on statistical 
questions revolving around the resolution of a spectrum estimate. It is 
generally known that, no matter what spectrum estimation technique is 
used, a decision must be made concerning the trade off between the resolu¬ 
tion of the estimate and its stability, or variance. The cost of increased 
stability is a loss in resolution, or vice versa. It is the possibility 
of taking advantage of the improved stability of the lower resolution 
spectrum estimate which makes it desirable to have a method of extracting 
a lower resolution estimate directly from a high resolution estimate. 

The ability to reduce resolution quickly and efficiently could make it 
possible to do frequency domain signal processing at different bandwidths 
without the expense of reprocessing the original time series. (In many 
cases the original time series is neither available nor recoverable for 
reprocessing.) An example of a situation in which this capability could 
prove useful is one in which signal spectra of various widths are to be 
detected in the presence of noise. In this case it would generally be 
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desirable to process the data originally with resolution comparable to 
the bandwidth of the narrowest line to be detected. Then, in order to 
enhance the detection capability for broader lines, the original high 
resolution spectrum could be altered to reduce its resolution and improve 
its stability. Using the techniques considered in this report, this 
could be done at much less computational expense than reprocessing the 
original time series. These techniques could also find application in 
the conversion of standard fixed resolution spectra to variable resolu¬ 
tion where, for example, resolution decreases with increasing frequency. 

The frequency averaging techniques to be examined here are not new, 
and the author is aware of several instances where they have been employed. 
In brief, they consist of the following procedure. If the resolution 
of the spectrum is to be reduced by a factor of four, for example, a 
number of high resolution spectrum values centered around the frequency 
of the desired low resolution value are averaged together. This is done 
for each low resolution spectral point desired. In the procedures to 
be examined here, the number of high resolution points combined to form 
a low resolution point is allowed to vary, and in some cases the averaging 
is done with unequal weighting of the high resolution values. 

In spite of their widespread usage, no detailed systematic examination 
of the simulated low resolution spectra produced by these techniques has 
been published. This report will make several comparisons between the 
simulated low resolution spectra and conventionally derived low resolution 
spectra. In particular, the questions to be addressed are: 
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(1) What are the statistical characteristics of a simulated low 
resolution spectrum estimate for a white noise process, and how do 
they compare with a conventional spectrum estimate which uses the 
same amount of data? 

(2) When the time series contains a narrowband signal embedded in 
white noise, how do the simulated and conventional spectra compare 
when both are derived from the same time series? 

(3) How do the simulated and conventional spectra differ in terms 
of minimum detectable signal level? 

(4) How does nonuniform weighting of the high resolution values 
in the frequency averaging affect the comparison in questions 1-3? 

What weighting functions yield superior performance? 

Question (1) is addressed analytically in Chapter III where expressions 
for the number of equivalent degrees of freedom for the simulated spectra 
are presented for the case of white Gaussian noise. These expressions 
are obtained for various direct spectrum estimation techniques including 
non-overlapped and 50% overlapped, averaged power spectra derived from 
discrete transformed time series, windowed with various data windows 
including the Hanning window. With some slight modification the techniques 
used in the derivations can be extended to other amounts of overlap and 
to other data windows. However, the cases examined here are representative 
of techniques currently being employed in practical situations. 

Questions (2) through (4) are not readily answered analytically 
since they are heavily data dependent. Instead, an empirical study via 
computer simulation is employed, and the results are presented in Chapter IV. 
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The results demonstrate that the weighting function used in averaging 
the high resolution values has a significant effect on the simulated 
low resolution spectrum. A procedure which yielded consistently good 
results is given for choosing the weighting coefficients. These 
coefficients produce simulated spectra whose response to narrow band 
signals is approximately the same as the conventional spectra, and which 
have approximately the same stability (as measured by equivalent degrees 
of freedom) as the conventional spectra. 

Of the remaining chapters. Chapter V contains a brief summary, 
while Chapter II is an introduction to, or a basic review of, spectrum 
estimation. The classic reference for this material is Blackman and 
Tukey [11 i. however, the text by Jenkins and Watts [ 2] gives a more 
modern and detailed treatment. 
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II. A BASIC REVIEW OF THE STATISTICAL NATURE \ 

OF ESTIMATES OF POWER SPECTRA j 


II.1 Discussion 

This chapter will serve as a brief review of some of the basic 
statistical facts associated with the estimation of power spectra from 
random data. It is not intended to be complete in any sense; however, 
it should orient the reader to the notation and viewpoint to be used 
in the remainder of the report. An excellent discussion of the entire 
subject is contained in Jenkins and Watts [2]. 

This report will be concerned with the performance of spectrum 
estimates when the available data x(t) consists of either a sample 
function taken from a stationary white Gaussian noise process, or such 
a sample function to which a sinusoid has been added. The stationary 
random process has an autocovariance function c x (t) defined by 

= E [x(t+x)-m x ] | , (2.1) 

where 

m x = Ejx(t)| . (2.2) 

The autocovariance function serves as a useful description of the process 
although it is not generally a complete description since it involves 
only second order moments. (It is a complete description, however, 
when the process is Gaussian.) A more widely used and somewhat more 
intuitively appealing characterization is the Fourier transform of the 
autocovariance function, or the power spectral density S x (f) of the 





process. It is the function S x (f) which we are interested in estimating, 
basing our estimate on an observation of a single realization of the 
random process over a time interval of length T. 

There are two basic approaches to forming estimates of S x (f). 
Until recent years, the more widely used technique consisted of using 
the observed x(t) to estimate c x (t) for discrete values of t, and then 
performing a discrete Fourier transform of the estimate of c x (t) to 
produce an estimate of S x (f). This is known as the indirect method, 
and a wide variety of techniques for implementing it have been discussed 
in the literature. More recently, with the advent of the FFT (fast 
Fourier transform) algorithm and faster computers, the preferred technique 
has become the direct method wherein x(t) itself is sampled and discrete 
transformed to the frequency domain where its average power, or mean 
square value, can be resolved into frequency components. The direct 
method is the one with which we will be concerned. 

As with the indirect method, the literature contains a wide 
range of techniques for implementing the direct method. With each of 
these techniques the objective is to use the available data to produce 
an estimate of the spectrum which has small bias, small variance, and 
the desired resolution. The bias of an estimate is its average error, 
the variance is a measure of the fluctuation of the estimate around 
its mean value, and the resolution of the estimate is determined by 
its ability to resolve spectrum components which are close together in 
frequency. Designing a spectrum estimator generally involves deciding 
on trade-offs between these three quantities, since it is usually 
possible to improve one at the expense of the others. We will be 
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primarily interested in the trade-off between variance and resolution. 

An example of how this trade-off occurs can be seen in a technique which 
is being used quite extensively to estimate spectra. In this technique 
the total data length T is broken up into shorter segments of length L 
which may or may not overlap each other (50% overlap is used commonly). 

The resulting segments are then sampled, transformed via the FFT, and 
a power spectrum estimate is obtained for each segment. These estimates 
are then averaged to form the final estimate. Variance reduction occurs 
because of the averaging; however, by segmenting the data and doing 
shorter transforms, we have decreased the resolution capability of the 
estimate. 

In the remaining sections of this chapter various examples of 
spectrum estimation will be discussed in somewhat more detail. A basic 
knowlege of probability and statistics is assumed. 

II.2 A Basic Estimate of the Spectrum 

Let the discrete time (sampled) version of x(t) be denoted 

x^ = x(kA) , 1 < k < L , (2.3) 

where 

A = l/f s (2.4) 

is the time interval between samples. Also, let 

L 

Aj = x^ cos(2TTjk/L) , (2.5a) 

k=l 

L 

B. = x k sin(2Trjk/L) , (2.5b) 

k=l 
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X. = A. + iB. , i = /-T 
3 3 3 


(2.5c) 


so that (X.) is the discrete Fourier transform of the sequence {x,}. 

3 x 

The sample spectrum based on x^, l<k^L, is generally defined as 


C (f.) = La)X. 
x j ' j 


( 2 . 6 ) 


where 


f. = (j-l)/LA 


(2.7) 


It can be shown [ 21 that under certain conditions C (f.) is an estimate 

x 

of S (f.) in the sense that 
x J 


lim e{c (f )} = S (f ) 

A-koo ' A ■> ' A J 


( 2 . 8 ) 


LA-*=° 


However C (f.) is not in itself a good estimate because it does not 
x J 

converge in any statistical sense to S (f.) as the record length L tends 

x 3 

to infinity. To illustrate this, we can examine the behavior of the 

sample spectrum for the case where the x^ are independent identically 

distributed zero mean Gaussian random variables (e.g., is a discrete 

white Gaussian noise process). Referring to Eq. (2.5) we see that and 

B. are linear combinations of Gaussian random variables, and so A. and 
J J 

Bj are themselves Gaussian. In fact, if is real Gaussian with mean 
2 2 

zero and variance o (i.e., x^ 'v, N(0,o )), then 


Aj, B. a, N(0, £-) , j = 1, \ 


(2.9a) 


A., B. a, N(0, , 1 < 3 < j 


(2.9b) 
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with the appropriate symmetry for values of j greater than j. So, for 
l<j<j, we have 


C x (fj) = LA(A?+B?) 



( 2 . 10 ) 


where A1 and Bl are N(0,1). The distribution of the sums of the squares 
of two N(0,1) random variables is well known to be the Chi-square dis¬ 
tribution with two degrees of freedom. Hence, for l<j<^ , 


C (f.) ^ 
x J 


2 . , 

o A 2 

2 X 2 


( 2 . 11 ) 


The mean of a Chi-square random variable with v degrees of freedom is 

simply v, while its variance is 2v, so for l<j<y , the mean of our 

estimate C (f.) is 
x 


E { C x (f j)> ■ ° 2t 

and 

Var|c x (f j )} = a 4 A 2 


( 2 . 12 ) 


(2.13) 


Note that the distribution of C (f.) does not depend on L, the amount 

* J 

of data processed. This result is not surprising when we realize that 
in this case where L is the transform length, increasing L increases 
our information about the spectrum, but it also increases the number 
of frequencies over which the information is spread. The result is no 
net gain in information at any one frequency, hence the variance remains 
constant. 
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II.3 Smoothing by Averaging 

Suppose now that we have a data record of length T sec and 
we wish to produce an estimate with variance somewhat smaller than that 
of the unsmoothed estimate in the preceeding section. This can be 
accomplished by segmenting the data into P separate pieces (no overlap), 
each of length LA. We can then form estimates C^(f-), l^p<P, in the 
manner described above, where each depends on a different segment 
of the data. The smoothed estimate is then the average, or 


p=l 


C?(f.) 


(2.14) 


The distribution of C is easily obtained from the fact that the sum of 

x 

independent Chi-square random variables has a Chi-square distribution 
whose number of degrees of freedom is the sum of the degrees of freedom 
of the component random variables. Hence 


C (f.) 
x j 


o 2 A 2 
2P X 2P 


(2.15) 


2 4 2 

with mean a A and variance a A /P. We have reduced the variance by a 
factor 1/P simply by limiting the resolution to P times what it would 
have been had we used a single transform. 

The analysis shown here is for the simplest type of smoothed 
spectral estimate. The use of overlapped transforms and data windows 
other than the simple rectangular one used here can result in greater 
variance reduction, but the statistical analysis can become complex. 
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II.4 The Use of Data Windows 

The data window implied in the above analysis is known as the 
rectangular window in which each L point FFT is performed on an unmodified 
segment of data. This is equivalent to multiplying the whole data record 
by a window function which is equal to unity for the segment to be 
transformed, and zero elsewhere. The resulting transform is then the 
convolution of the desired transform with the transform of the window 
function, which in this case is a sine function (with a phase factor). 
Therein lies the chief drawback of the rectangular data window. The 
sidelobes of the sine function are large, implying that, for spectrum 
components whose period is not an integral divisor of LA, there can be 
significant "leakage" or spreading of energy into adjacent frequencies. 

The leakage problem can be reduced by using another data window whose 
transform has smaller sidelobes, such as the Hanning or cosine window 
(see Appendix 1). The window can be applied in the time domain by 
multiplying each data value by the appropriate window value before 
performing the FFT. One advantage of the Hanning window, however, is 
that it can also be applied quite efficiently in the frequency domain 
via a three point convolution. 

It is interesting to consider the effect of the Hanning window 
on the variance of an averaged spectrum estimate such as the one discussed 
in Section II.3. As is pointed out in Appendix I, the window can be 
applied by convolving the complex coefficients resulting from the FFT 
of the rectangular windowed time sequence with the sequence 
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W(k) 


0.5 , k = 0 

-0.25 , k = ±1 

0 , |k| > 1 


(2.16) 


If we allow A^ ^ and ^ to denote the real and imaginary parts of the 
Hanned complex coefficients, and let A^ and be as defined in Eq. (2.5), 
then 


A, . * -0.25 A. , + 0.5 A. -0.25 A. , (2.17a) 

n,J j-1 3 J+l 

B. . = -0.25 B. . + 0.S B. -0.25 B. , . (2.17b) 

h,J J-1 j J+l 

Again allowing the sequence to be a realization of a discrete 
white Gaussian noise process, we can easily show that for 2<j<^-2, 

2 

Aj^ -w N(0, 0.375 |jj) , (2.18a) 

2 

B h . i N(0, 0.375 |j-) . (2.18b) 


Then the distribution of the spectral estimate obtained by averaging P 
independent Hanned transforms can be obtained using the same techniques 
used previously: 


C (f.) * 0.375 
x j 


2 , 
o A 2 

2P X 2P 


(2.19) 


We see that for the white Gaussian noise case, (f^) is essentially 

the same random variable as C (f.), multiplied by a constant scale factor. 

x J 

(Alternatively, we could note that C x (f^) and C x (f^) have the same number 
of equivalent degrees of freedom.) So, although the Hanning window does 
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reduce the leakage problem, it actually does not contribute to estimate 

stability in the case of averaged non-overlapped transforms. 

As a final comment on the use of the Hanning window, we would 

like to point out its effect on the frequency autocovariance of the 

spectrum estimate. The autocovariance function is derived in Appendix 3 

and will only be stated here. Let be a spectrum estimate resulting 

N 

from a single Hanned transform. Then, for 2<i,j<j -2, 


Cov{c-Cfi), W} A 


1.406*10" Va 4 , | i-j | = 0 


6.250xl0' 2 o 4 A 2 , ji-j | = 1 


3.907xl0' 3 a 4 A 2 , |i-j| =2 


( 2 . 20 ) 


0 . 

V 


-j I * 3 


For P non-overlapped Hanned transforms. 


Cov 


{ c x (f i )> - F 


'fc' 


c ° v ( c ;< f i>- c ;<v} 


( 2 . 21 ) 


From Eq. (2.21) we see that averaging several non-overlapped Hanned spectra 
has no real effect on frequency autocovariance since normalization of the 
resulting functions yields the same normalized autocovariance for both. For 
comparison purposes we note that when the rectangular data window is used 
in the white Gaussian noise case, the Fourier coefficients are independent 
resulting in 


4 a 2 
a A 


,i = J 


c ° v { c x <f i>' C x <f j>H 


0 , i ? j 


( 2 . 22 ) 







We can conclude that although leakage is reduced, which is important 
when narrow band signals are present in the data, we have introduced 
some correlation in the noise portion of the spectrum by using a window 
other than the rectangular window. 

II.5 Equivalent Degrees of Freedom 

There are a number of spectral estimation techniques which 
yield estimates whose distributions, unlike those discussed above, are 
not that of a constant times a Chi-square random variable. To compare 
the stability of these estimates, a common approach is to approximate 
the perhaps unknown distribution of the spectrum estimate with a 
distribution which i£ that of a constant times a Chi-square random 
variable. The approximation is in terms of first and second order 
moments. Let X be the random variable whose distribution is to be 
approximated by that of aY, where Y is Chi-square with v degrees 
of freedom. We assume that the mean and variance of X can be calculated. 
Then, in order to determine appropriate values for a and v, we equate 
means and variances: 


(2.23) 

(2.24) 

Solving for a and v yields 


EX = EaY = otEY = av 


Var X = Var aY = a Var Y 


, 2 

= 2a v 


Var X 

01 = ~2Er 


2 

2E X 
Var X 


(2.25) 

(2.26) 
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Hence, if we can compute the mean and variance of the spectrum estimate, 
we can compute the number of degrees of freedom of a Chi-square approxima¬ 
tion to its distribution. This is termed the number of equivalent degrees 
of freedom of the estimate. Its value as an indicator of the stability 
of an estimate lies in the fact that it increases as the size of the 
variance, relative to the mean, decreases. It should be borne in mind, 
however, that the Chi-square distribution corresponding to the equivalent 
degrees of freedom is only an approximation to the true distribution of 
the spectrum estimate. 

II.6 Summary 

In this chapter and the appendices referred to, we have tried 
to present some basic material which is fundamental to the remainder of 
the report. Both the material presented and the presentation are by no 
means complete; however the interested reader should be able to satisfy 
the shortcomings by referring to Jenkins and Watts [2 ], Bergland [3], 
Cooley et al. [41, or Richards [5]. These references are only a 
sampling of a number of papers and texts addressing basic topics related 
to direct spectrum estimation via the FFT. 




III. STATISTICAL STABILITY OF LOW RESOLUTION SPECTRA 
OBTAINED FROM FREQUENCY AVERAGED HIGH RESOLUTION SPECTRA 


III.l Discussion 

The frequency averaging methods by which low resolution spectrum 
estimates are derived from high resolution estimates are straightforward 
and intuitively obvious. We simply partition the frequency range into a 
number of intervals or bands whose widths correspond to the desired fre¬ 
quency resolution for the low resolution estimate. We then perform a 
weighted average over all the high resolution spectrum values in each 
band. The procedure can be generalized by allowing the low resolution 
bands to overlap so that some of the high resolution points may enter 
into more than one of the weighted averages. 

In this chapter we will present formulae which will allow 
computation of the number of e_quivalent degrees of freedom of a spectrum 
obtained by frequency averaging. We will assume that the time series 
on which the estimates are based is simply a discrete white Gaussian 
noise process. The derivations have been completed for cases where the 
high resolution spectrum was estimated using both overlapped and non- 
overlapped transforms, with rectangular, Hanning, and other data windows. 

The principal importance of these results is that they allow us to 
determine weighting coefficients for frequency averaging which will 
produce a low resolution estimate with the desired stability. For example, 
in some of the empirical studies to be presented in Chapter IV, the weighting 
coefficients were chosen to produce a simulated low resolution spectrum 
with the same number of equivalent degrees of freedom as a conventionally 
obtained spectrum based on the same amount of data. 
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The results are grouped in Sections III.2-III.4 according to 


the method used to produce the high resolution estimate. The derivations 
are contained in Appendices 2 and 4. 

III.2 Non-overlapped High Resolution Transforms with the Rectangular 
Data Window 

Let C(fj) denote a power spectrum estimate produced by averaging 
P non-overlapped rectangular windowed transforms, and define the frequency 
averaged low resolution spectrum derived from C by 


e 'V ■ £ “j 

j = -K 


<C( f k .;j> 


(3.1) 


The low resolution spectrum estimate for frequency fj is then a weighted 
sum of 2K+1 high resolution spectrum values centered around f . The 
computation of the number of equivalent degrees of freedom v of C(f k ) 
can be carried out easily using Eq. (2.26). 


v = 


2E 2 C(f k ) 
Var C(f R ) 


(3.2) 


2E2 i E V (f k-J 3 


>J = ~K 


a j C(f k-j : 

J = ~K 


It was pointed out in Chapter II that the random variables C(f^) are 
independent over j, so we have 


(3.3) 
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V 


(3.4) 



2 a ' var f''k-i 1 } 

j = -K 


Also from Chapter II, we know that for values of j sufficiently removed 
from the endpoints of the spectrum. 



(3.5) 

(3.6) 

(3.7) 


III.3 Non-overlapped High Resolution Transforms with the Hanning 
Data Window 

For this case, only the results will be stated and the derivation 


will be deferred to Appendix 2. Using the same notation used above. 
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where 


D(k.j) 


^0.4219 , 

k=j 


0.26S6 , 

1 k-j | 

= 1 

0.1719 , 

Ml 

= 2 

0.1406 , 

V. 

Ml 

Al 


f 3.9) 


In cases where the weighting coefficients are all equal and the 
frequency averaging includes N>4 high resolution spectrum points, 
Eq. (3.8) reduces to 


S.63N 2 P 
S.47N - 2.81 


(3.10) 


III.4 50% Overlapped High Resolution Transforms 

For this case also, only the results will be stated. The 
derivation is in Appendix 4. In this case. 



( 3 . 11 ) 


( 3 . 12 ) 

(3.13) 
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(4) 

(5) 

( 6 ) 



(3.15) 


P is the number of overlapped high resolution transforms 
averaged together to form C Cfj D > and 

W(f) is the spectral window, or the Fourier transform of the 
data window. 

When the Hanning window is used (Appendix 1), 


W(f) = L sinc ( Lf ) 
2[l-(Lf) 2 ] 


(3.16) 


Evaluating the integral in Eq. (3.12) yields, for the Hanning window, 

I Q = 0.1406L 2 . (3.17) 

The integrals in Eqs. (3.13) and (3.14) have been evaluated numerically 
for the Hanning window: 


f 

1.406 x 

10 _1 L 2 , 

k=j 


6.250 x 

io' 2 l 2 , 

|k-j| 

= 1 

3.906 x 

io' 3 l 2 , 

|k-j 1 

= 2 

= 0 

} 

1 k-j | 

IV 

c Si 


V. 


(3.18) 
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3.906 x 10-V 


2.814 x 10-V 


|l 2 (k.j)|" = ^ 9.766 x 10'V 

-4 2 

1.126 x 10 L 


« 0 


k=j 

|k-j| = 1 

Ik-j| = 2 

I k-j| = 3 
|k-j| * 4 


Substituting these values back into Eq. (3.11), we have for the 
Hanning window, 


v = 


2.813 x io 1 p| ^ a_. 

' j = -K 


H, 


k=0 


£ Vj 

i» j : 

i-jl=k 


where 


H = 1.406 x 10 _1 + (3.906 x 10' 3 ) 

O r 


H 1 = 6.250 x io' 2 + (2.814 x 10‘ 3 ) 


H 2 = 3.906 x 10 


-3 + 2(P-1) 


(9.766 x 10 ) 


H 3 = 2(P p" 1) f 1 - 126 x 10 " 4 ) 
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(3.19) 


(3.20) 


(3.21a) 
(3.21b) 
(3.21c) 
(3.21d) 






Ill.5 Examples 


In this section we will use the above results to calculate 
equivalent degrees of freedom for frequency averaged low resolution 
spectrum estimates. The calculations will be done for two types of 
weighting coefficients used in the averaging. Since the most commonly 
used coefficients are uniform (all the coefficients are equal), this 
case was chosen for one example. In addition, some effort was devoted 
to finding a method of choosing coefficients which would produce a 
frequency averaged low resolution spectrum whose characteristics were 
as similar as possible to those of a conventionally obtained spectrum 
with the same resolution based on the same data. This resulted in a 
set of weighting coefficients which will be referred to as the HSQ 
coefficients, obtained from the following equation: 


a. 
J 


sine (aj) 

f . 2 . 2,2 

(1 -“) ) 


(3.22) 


Here, a is the ratio of the desired low resolution to the higher 
resolution from which the frequency averaged spectrum is to be obtained. 
The coefficients derived from Eq. (3.22) can then be used in Eq. (3.1) for 
the frequency averaging. 

While doing the examples for this report, we were primarily 
interested in reducing resolution by factors of two, four, and eight, 
corresponding to a values of 0.5, 0.25, and 0.125, respectively. The 
example computations are performed for a situation where the high 
resolution transform length L was one-sixth of the total available 
data T. Hence P was equal to six if non-overlapped transforms were 
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used, or eleven if 50% overlapped transforms were used. The 
results of the calculations are shown in Figs. 1-4 for the four types 
of weighting coefficients. In each figure are shown three curves 
corresponding to high resolution processing using non-overlapped trans¬ 
forms with rectangular and Hanning windows, and 50% overlapped 
transforms with the Hanning window. The number of equivalent degrees 
of freedom of the frequency averaged spectrum was plotted versus the 
number of coefficients used. For example, if K were equal to seven in 
Eq. (3.1), then 15 coefficients would be used. If these coefficients 
were uniform, and if the high resolution spectrum was derived from non- 
overlapped Hanned transforms, then, from Fig. 1, the low resolution 
spectrum would have approximately 95 equivalent degrees of freedom. 

The effect of the Hanning window in introducing correlation 
across frequency into the spectrum appears in each example. This is 
the reason why the rectangular windowed spectrum, with its absence of 
frequency correlation, produces more equivalent degrees of freedom per 
window coefficient than the Hanning windowed spectrum. This is 
generally true except when overlapped transforms and a small number of 
coefficients are used. In this case the effect of the increased number 
of transforms outweighs the effect of frequency correlation resulting 
in a larger number of equivalent degrees of freedom than produced by 
non-overlapped transforms with the rectangular window. 

For comparison purposes, we point out that from [Ref. 6, 

Eq. (36)J we can compute the number of equivalent degrees of freedom for 
the high resolution spectrum when 50% overlapped Hanned transforms are 
used. The resulting value is approximately 21. A 2:1 reduction 
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EQUIVALENT DEGREES OF FREEDOM 



FIGURE 2 

EQUIVALENT DEGREES OF FREEDOM versus WINDOW LENGTH 
FOR 2:1 HSQ COEFFICIENTS, T - 6L 
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FIGURE 3 

EQUIVALENT DEGREES OF FREEDOM versus WINDOW LENGTH 
FOR 4:1 HSQ COEFFICIENTS, T = 6L 
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FIGURE 4 

EQUIVALENT DEGREES OF FREEDOM versus WINDOW LENGTH 
FOR 8:1 HSQ COEFFICIENTS, T = 6L 


ARL:UT 
AS-78-1490 
CSP • GA 

fi.1K.B1 








(*• 


in resolution with the same type of processing using shorter transforms 
and the same amount of data would yield 44 equivalent degrees of freedom, 
while a 4:1 reduction would yield 89 and an 8:1 reduction would yield 
180 equivalent degrees of freedom. For the examples to be presented in 
the following chapter we will choose a window length which produces 
approximately the same number of equivalent degrees of freedom as would 
have been produced by conventional low resolution processing of the same 
data. In other words, if we were to take the high resolution spectrum 
in this example and reduce the resolution by frequency averaging instead 
of reprocessing the data, we would choose a window length based on the 
graphs of Figs. 1-4 which would yield 44, 89, or 180 equivalent degrees 
of freedom for 2:1, 4:1, or 8:1 reduction ratios. This would insure 
that the frequency averaged spectrum would have the same statistical 
stability as a spectrum produced by reprocessing the original data with 
shorter transforms. 


It 
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CHAPTER IV 


EFFECTS OF FREQUENCY AVERAGING ON NARROWBAND SIGNALS 

IV.1 Discussion 

In the preceeding chapters, interest has been primarily in 
the estimation of power spectra for data consisting solely of random 
noise. In Chapter III, formulae were obtained which would allow 
determination of the number of equivalent degrees of freedom possessed 
by a frequency averaged low resolution spectrum as a function of the 
averaging coefficients. In this chapter, we will be concerned with 
the performance of the averaging techniques when the data contains a 
sinusoid in addition to the random noise. 

Our objective here is to compare various aspects of the 
performance of spectrum estimates obtained by frequency averaging with 
estimates obtained by conventional processing of the data, where the 
comparisons are made with effects on detection in mind. The comparisons 
are made by generating a data sequence consisting of a sinusoid plus 
random noise. Two spectrum estimates are then produced from the same 
data. One estimate, the conventional one, employs fifty percent over¬ 
lapped FFTs in the manner described previously. The other estimate, 
the frequency averaged one, is derived from a higher resolution 
spectrum which also resulted from fifty percent overlapped FFTs. In 
all of the examples presented here, the weighting coefficients are 
chosen to produce a frequency averaged spectrum with approximately the 
same number of equivalent degrees of freedom as the conventionally 
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derived spectrum. Examples are presented in which the resolution 
reduction is by factors of 2:1, 4:1, and 8:1. Spectrum comparisons 
are made for several sinusoidal signal frequencies evenly spaced 
between FFT frequencies. The weighting coefficients used are the HSQ 
and uniform coefficients described in Section III.5. The coefficients 
used in all the examples in this chapter are shown in Table 4.1. 

They are normalized so that their sum is unity. 

IV.2 Comparison of Spectral Windows 

This section is concerned with the shape of the spectral 
window associated with frequency averaged spectra. Knowlege of the 
spectral window, and how it compares with such commonly used windows 
as the Hanning window, is important since it indicates how much leakage 
of energy into adjacent frequencies will occur for sinusoids whose 
period is not an integral divisor of the time length of the FFT. 

Spectral windows were determined for signals in the absence of noise 
for both HSQ and uniform coefficients. They are plotted in Figs. 5-7 
along with the Hanning window used in producing the conventional estimate. 

A few words of explanation are in order concerning these 
figures. The abscissa values of frequency are normalized by multiplying 
the actual signal frequency by the time length of the low resolution 
transform used in the conventional estimate. Hence integer values of 
frequency correspond to the discrete frequencies at which the FFT 
evaluates the transform. Spreading can be determined from the graphs 
as follows. Suppose we have a signal which is 0.75 normalized frequency 
units away from an FFT frequency. If the Hanning window is used with 
conventional processing, it will produce a contribution to the energy 
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TABLE 4.1 

COEFFICIENTS FOR FREQUENCY AVERAGING 


2:1 Reduction 


HSQ Coefficients 


K = 2, EDF = 46.5 

a = 0.3400 
o 

a +1 = 0.2450 
a +2 = 0.0850 


Uniform Coefficients 


K = 1, EDF = 38.4 

a = 0.3333 
o 

a , - 0.3333 
+ 1 




K = 6, EDF = 89.3 

a = 0.1670 
o 

a = 0.1540 

a +2 = 0.1203 

a +3 = 0.0786 

a +4 = 0.0417 

a +s = 0.0171 

a,, = 0.0048 
±6 


K = 12, EDF = 174.6 


K = 4, EDF = 100.5 

a =0.1111 
o 

a = 0.1111 

a +2 = 0.1111 

a +3 = 0.1111 

a . = 0.1111 
±4 


K = 8, EDF = 184.3 


a = 0.0836 

a = 0.0588 

0 

0 

a , = 0.0819 

a , = 0.0588 

+ 1 

±1 

a „ = 0.0771 

a _ = 0.0588 

±2 

±2 

a , = 0.0696 

a „ = 0.0588 

+3 

±3 

a ± . = 0.0602 

a . * 0.0588 

±4 

±4 

a iC = 0.0498 

a c ' 0.0588 

+5 

±5 

a,, = 0.0393 

a , = 0.0588 

±6 

±6 

a^_ = 0.0295 

a _ = 0.0588 

±7 

±7 

a „ = 0.0209 

a 0 = 0.0588 

±8 

±8 

a „ = 0.0139 


±9 


: = 0.0086 


±10 


. .. = 0.0048 


±11 


. , = 0.0024 


±12 

































FIGURE 7 

SPECTRAL WINDOW COMPARISONS, 8:1 REDUCTION 
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at that FFT frequency which is 3.4 dB less than its contribution would 
have been had its frequency been the same as the FFT frequency. 

The graphs reveal that HSQ coefficients, matched to conventional 
processing in equivalent degrees of freedom, also produce a spectral 
window which matches the Hanning window quite closely. The graphs also 
reveal that when enough uniform coefficients are used to produce the 
required number of equivalent degrees of freedom, the resulting spectral 
window is somewhat broader resulting in less resolution capability. 
Resolution can be restored by using fewer coefficients, but this will 
be at the cost of spectrum stability in terms of equivalent degrees of 
freedom. Finally, the graphs show that the spectral window can be 
tailored significantly by varying the coefficients used in the frequency 
averaging. 

IV.3 Comparison of Estimated Spectra 

Figures 8-25 are plots of various estimates of the power 
spectrum of a sinusoid in the presence of noise. A few words of 
explanation are in order concerning these plots. A 1 Hz sampling rate 
was used, so the spectra are shown from DC to the Nyquist frequency of 
0.5 Hz. The sinusoid (signal) frequency is given in the title of each 
plot in normalized frequency units, i.e., the true frequency is multiplied 
by the length of a low resolution (short) transform. In each of these 
examples, the low resolution transform length was chosen to be 128 points, 
while the high resolution transforms were 256, 512, or 1024 points for 
2:1, 4:1. and 8:1 reductions in resolution. Using the normalized fre¬ 
quency units, integer values of signal frequency correspond to FFT 
frequencies. Since the response of the spectrum estimates to sinusoids 
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2:1 REDUCTION, L x F * 32.000 


ARL:UT 
AS-78-1494 
CSP - GA 
6-15-81 


38 






HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 
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FIGURE 9 

COMPARISON OF POWER SPECTRA 
2:1 REDUCTION, LxF* 32.100 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


FIGURE 10 

COMPARISON OF POWER SPECTRA 
2:1 REDUCTION, L x F = 32.200 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 
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FIGURE 11 

COMPARISON OF POWER SPECTRA 
2:1 REDUCTION, LxF - 32.300 













HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


FIGURE 12 

COMPARISON OF POWER SPECTRA 
2:1 REDUCTION, L x F = 32.400 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


FIGURE 13 

COMPARISON OF POWER SPECTRA 
2:1 REDUCTION, LxF = 32.500 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


FIGURE 14 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, L x F = 32.000 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 
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FIGURE 15 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, Lx F = 32.100 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 
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FIGURE 16 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, L x F = 32.200 
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FIGURE 17 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, L x F = 32.300 
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APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 
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FIGURE 18 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, L x F = 32.400 









HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 
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FIGURE 19 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, L x F = 32.500 
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HIGH RESOLUTION SPECTRUM 
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FIGURE 20 

COMPARISON OF POWER SPECTRA 
8:1 REDUCTION, L x F = 32.000 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


ARL:UT 
AS-78-1508 
CSP - GA 
6-15-81 


FIGURE 22 

COMPARISON OF POWER SPECTRA 
8:1 REDUCTION, Lx F = 32.200 
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FIGURE 23 

COMPARISON OF POWER SPECTRA 
8:1 REDUCTION, LxF = 32.300 
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HIGH RESOLUTION SPECTRUM 
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FIGURE 24 

COMPARISON OF POWER SPECTRA 
8:1 REDUCTION, L x F = 32.400 
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HIGH RESOLUTION SPECTRUM 
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FIGURE 25 

COMPARISON OF POWER SPECTRA 
8:1 REDUCTION, LxF = 32.500 
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will be symmetric in frequency around each FFT frequency, a representa¬ 
tive idea of the effect of signal frequency on the behavior of the 
spectrum estimates can be had by varying normalized signal frequency 
between 32.0 and 32.5, for example. This was done in steps of 0.1 
normalized frequency units for 2:1, 4:1, and 8:1 reduction ratios in 
the examples. 

The amplitude normalization for the examples was chosen to 
preserve the amplitude of the signal for both the high and low resolu¬ 
tion estimates. Since the noise power is uniformly distributed over 
frequency while the signal power is concentrated at a single frequency, 
increasing resolution spreads the noise power over a larger number of 
FFT frequencies. The signal, on the other hand, remains concentrated 
in a small number of FFT frequencies so the net result is a decrease 
in the amplitude of the mean noise level relative to the signal level. 
Since the spectra are normalized to preserve signal levels in the examples, 
the mean noise level in the low resolution estimates will be higher than 
in the high resolution estimates. 

In each figure, the top plot is the high resolution spectrum 
from which the frequency averaged spectra are derived for both HSQ and 
uniform coefficients. In the middle plot are shown both the low resolu¬ 
tion spectrum obtained by reprocessing the data with shorter transforms, 
and the frequency averaged spectrum obtained with HSQ coefficients. In 
the bottom plot, the same low resolution conventional spectrum is shown 
with the frequency averaged spectrum obtained using uniform coefficients. 
From the figures, it is quite clear that both sets of coefficients produce 
spectra which compare quite closely with the conventionally obtained 
spectrum. 
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IV. 4 


Signal Excess Comparison 

To get a first order approximation of the gain or loss in 
ability to detect narrowband signals when frequency averaged spectra 
are used in place of conventionally obtained spectra, signal excess 
values were calculated for the examples. These values are 
presented in Figs. 26-28. They were obtained in the following manner. 

First, the high resolution spectrum, the low resolution spectrum, and 
both HSQ and uniform coefficient frequency averaged spectra were computed 
for a data sequence. Signal excess in dB was determined for all three 
of the low resolution spectra by referring the peak value of the spectrum 
to the known mean of the noise portion of the spectrum. This was repeated 
ten times on independent data sequences at each of 40 signal frequencies 
between FFT frequencies. The resulting average values of signal excess 
as a function of signal frequency are then shown in the top plots of 
Figs. 26-28 for 2:1, 4:1, and 8:1 resolution reduction ratios, respectively. 
In the bottom plots of these figures are shown the difference in signal 
excess between the conventional low resolution spectrum and the frequency 
averaged spectra. 

The top plot in these figures shows the expected gradual decline 
in signal excess as the signal frequency approaches the midpoint between 
FFT frequencies. This is what is commonly referred to as the scalloping 
effect. In the bottom plot of Fig. 26 we see that the difference 
between the frequency averaged spectra and the conventional spectrum 
has the opposite shape, with a peak in the middle. This can be 
explained by examining the spectral windows in Fig. 5. Here it is seen 
that both of the frequency averaged spectral windows are somewhat 
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FIGURE 27 
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FIGURE 28 
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greater than the Hanning window at frequencies near the midpoint between 
FFT frequencies. We see that as the HSQ coefficient spectral window 
becomes a better approximation to the Hanning window at 4:1 and 8:1 
reduction ratios, the signal excess difference in Figs. 27 and 28 becomes 
quite small. On the other hand, the differences between the Hanning 
window and the uniform coefficient window seen in Figs. 6 and 7 result in 
obvious signal excess differences. This suggests that the performance 
of a set of coefficients to be used in frequency averaging can be tested 
quickly by determining the associated spectral window and comparing it 
with the desired spectral window. 
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IV.5 Summary 

In this chapter we have shown examples of frequency averaged 
spectra where our primary goal was to produce a lower resolution spectrum 
with characteristics as similar as possible to those of a low resolution 
spectrum obtained by reprocessing the original time series data. Both 
sets of coefficients yielded good results, with the HSQ coefficients 
holding a substantial advantage. In fact, the accuracy of the HSQ co¬ 
efficient estimate can be improved still further by adjusting the value 
of a in Eq. (3.22). 

We have pointed out the fact that the frequency averaging 
coefficients have a substantial impact on the effective spectral window 
of the frequency averaged spectrum. In Appendix 5 we show further 
examples where resolution reduction of 4:1 was done using four uniform 
coefficients. These examples, in which neither the number of equivalent 
degrees of freedom nor the effective spectral window match those of 
conventional processing, is included because it illustrates the most 
commonly employed technique. Some of the shortcomings of this technique 
are discussed in the appendix. 
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V. CONCLUSIONS 


The results and examples presented here are by no means the 
last word on frequency averaging to achieve resolution reduction. 

However, they should form a sound basis of understanding for both 
current usage of frequency averaging techniques and for further develop¬ 
ment and exploration of these techniques. 

As a basis for current usage, the report presents equations 
for calculating the number of equivalent degrees of freedom of frequency 
averaged spectra as a function of the coefficients used in the averaging 
process. Some guidance in the selection of coefficients is provided by 
the examples in Chapter IV and Appendix 5. Computer programs to numerically 
determine the spectral window associated with a set of coefficients can 
be easily written, enabling the use of the coefficients to tailor the 
window to suit particular needs. This capability for shaping the spectral 
window is itself a feature of the frequency averaging technique whose 
utility should be explored. Future work in this area should also examine 
the impact of signal bandwidth on these techniques. 
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APPENDIX 1 

NOTES ON THE HANNING DATA WINDOW 






Let 


rect(x) 


I< x < I 

2 - x _ 2 


1*1 > f 


Then the Hanning data window for a record of length L seconds extending 
from t= - j to t= j is defined as 


w(t) = j|l + cos rect(i] 


Let W(f) be the Hanning spectral window defined by 

W(f) = F(w(t)} , 

where F denotes the Fourier transform. 

Then, using the asterisk to designate convolution. 


W(f) = j F{1 + cos(2-rrt/L)} * F(rect(t/L)} 


H S(f) 


L sin(wLf) , L “"Wf-C» . 1 


wLf 


2irL(f-i) 


2nL(f+ j- ) 


sm 

2 


w) ri i f _L_ + l j 

it [f ■ 2 l Lf-l Lf+l ; J 


L sinc(Lf) 
2[l-(Lf) 2 ] 
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where sinc(x)=sin(ux)/nx. The usual procedure is to normalize the data 
window so that 

w ( f) = sinc(Lf - ) ? . 

1 - (Lf) 

In many cases of practical interest, the data window extends 
from t=0 to t=L: 

w’ (t) = w(t- y) 

Then W'(f) = F{w* (t)} 


-imLf 

= e 


W(f) 


Since discrete transforms are generally used to process data, we will be 
interested in W'(f) evaluated at f^=k/L for integer values of k: 


w- ( f k ) 


= < 


L 

2 

L 

'4 


k=0 

k= ± 1 
I k I > 1 


The convenient discrete form of W'(f^) allows us to implement the data 
window by convolving the rectangular windowed data with the sequence 
(- y, y, - y), or, where normalization is used, with (- y, y, - y). For 
example, if i=/-T and 


X. 

J 


A. + iB. 
1 1 


is the complex discrete transform of the sequence x^, where the rectangular 

k 

data window has been used, we can form the complex Hanned transform X., by 
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This is frequently more convenient than multiplying the data sequence 
by the discrete version of w(t) before transforming. 
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APPENDIX 2 


STATISTICAL STABILITY OF LOW RESOLUTION SPECTRA DERIVED FROM 
NON-OVERLAPPED HIGHER RESOLUTION FFTs 
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In this appendix we will assume that we have data which consists 
of a sample function from a discrete stationary zero mean white Gaussian noise 
process. From this data we will extract a number P of non-overlapping 
segments, each of which will be discrete transformed with the Hanning 
window applied. An estimate of the power spectrum of the process will 
be obtained from each discrete transform and the resulting P spectrum 
estimates will be averaged to form what will be referred to as the high 
resolution spectrum estimate. From this estimate, a lower resolution 
estimate will be obtained by means of the frequency averaging techniques. 

Our objective here is to obtain an expression for the number of equivalent 
degrees of freedom of the lower resolution estimate. 

Let C p(fj) be high resolution estimate obtained from the p tfl 
segment of data, where fj = (j-l)/LA, L being the length of the discrete 
transform and A being the sampling interval of the process. Then, the 
lower resolution estimate obtained by frequency averaging is 


C(fj) 



(A-2.11 


(A-2.2) 
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Since the random variables C^Cf^) are independent and identically 
distributed over lsp^P, the following lemma will prove useful- 
Lemma . Let Z^, .... Z p be independent identically distributed random 
variables, not necessarily Chi-square. Let 


2 

2E Z 


v = 


1 


Var Z 


1 


be the number of equivalent degrees of freedom of each Z^. Further, let 


■fEh 


i = l 

Then the number of equivalent degrees of freedom of Z is Pv. 
Proof. The number of equivalent degrees of freedom of Z is 


2E‘ 


2 

2E Z 

Var Z / P 

Var < ^ 1 

i = l 


I'S-i 

l-s-4 


2 

2E Z, 


P Var Z 1 


= Pv 


which completes the proof. 


74 





r 


We can now compute the number of equivalent degrees of freedom of 
C(fj) by computing the number of equivalent degrees of freedom for the 
bracketed term in Eq. (A-2.2) for p=l, and invoking the lemma. For 
simplicity, let 

Cffj) = CjCf.) (A-2.3) 


be the spectrum estimate based on the first segment of the data, or 

x., x . Let X., l<j<L, be the discrete transform of the data, 

l L j 

and 1 et 

X. = A. + iB. , (A-2.4) 

3 3 1 

where i = /^T. As pointed out in Chapter 2, the A. and B., l<j<L, are 

2 

independent zero mean Gaussian with a variance which we will denote a . 
Applying the Hanning window via a convolution in the frequency domain, 
we have 


1 


i. , + 0.5 A. 
J" 1 3 

- 0.25 A. . 

3+1 

CA-2.5a) 

i. . + 0.5 B. 
3*1 3 

- 0.25 B. , 

3 + 1 

(A-2.5b) 

3 


(A-2.5C) 


Then the estimate C(f^.) is given by 


C(f.) = LA|X. 


= LA(A 2 + 8 2 ) 
J J 
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(A-2.6) 








Without loss of generality, we can assume LA=1, since multiplication 
of a random variable by a constant does not alter its number of 


equivalent degrees of freedom. Also, by assuming that A , ..., A^ + j 

2 

and B q , .... B^ +1 are independent zero mean Gaussian with variance o , 
we can restate our objective as the calculation of the equivalent 
degrees of freedom of 



(A-2.7) 



* k = l ' 

-P- • (A-2.10) 

Var j£ a k A k I 

1 k=l ' 
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2 2 

Since J^a^A^ and ^a^B^ are independent and identically distributed. 
Expanding Eq. (A-2.10) yields 


v 


N N 


EE 

Lr - 1 n = l 




(A-2.11) 


since the A k are identically distributed. Some fairly tedious computation 
can be employed to show that 


EA. 2 = 3. 
k 

.750 x 10' 

1 2 
a 

1 s k < N 

(A-2.12) 


r 

-1 4 

il 



4.219 x 

10 a 



2.656 x 

O 

i—i 

l 

O 

, |k-j| = 1 


ea:a = < 



(A-2.13) 

k J 

X 

00 

iH 

10 _1 o 4 

CM 

II 

1 



1.406 x 

D 

f—i 

i 

O 

i-H 

, |k-j| > 3 





As an example, when uniform weighting c<' lci<- ‘ are used, 

{ 2 , N = 1 

2.77 , N = 2 , (A-2.14) 

3.72 , N = 3 


and, for N>4 in this case, Eq. (A-2.11) reduces to 

























APPENDIX 3 


THE EFFECTS OF THE HANNING WINDOW ON FREQUENCY AUTOCOVARIANCE 





Let be independent zero mean Gaussian random variables 

2 

with variance a , and let X^, .... be their discrete Fourier transform 
where 

X. = A. + iB. . (A-3.1) 

J J 3 ; 

As pointed out in Chapter 2, the random variables A. and B. are independent 

2 J J 

zero mean Gaussian with variance for . The corresponding spectrum 

estimate is 


c(f.) = la|x. r 

= LA(A? + B?) . (A-3.2) 

The frequency autocovariance is then 

Cov {C (f _.) , C(f.)} = E|[C(f.) - E{C(f j )}][C(f i ) - E{C(f.)}]| 

i 

i 

= E{C(f.) C(f.)> - E(C(f.)} E(C(f i )} . (A-3.3) 


Note that for i/j, C(f.) and C(f_.) 
is zero. If i=j. 


are independent, hence their covariance 


Cov{C(f ), C(f.)} = Var{C(f.)} 
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t- 




from Eq. (2.13). Hence, for a single transform estimate of the spectrum. 


Cov(C(f.), C(f.)) 



i = j 
i * j 


(A-3.4) 


In similar fashion, if C(f^) is the estimate which results from averaging 
P non-overlapped transforms, we can show 


Cov{C(f ), C(f.)) 



i = j 
i / j 


(A-3.5) 


Suppose now that we apply the Hanning data window via convolution 
(see Appendix 1). Then the complex Hanned transform is 

X. = A. + iB. , (A-3.6) 

J J J 

where 


A. - 

-0.25A. , 

+ 0.5A. 

- 0.25A. . 


(A-3,7a) 

J 

J-l 

J 

J+l 

9 

8. = 

J 

-0.25B. , 

J-l 

+ 0.5B. 

J 

- 0.25B. . 
1 + 1 


(A-3.7b) 


The Hanned spectrum estimate is 

C(f.) = LA | X j 1 2 

= LA(A 2 + 8 2 ) . (A-3.8) 

j y 
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As before, the frequency autocovariance is 


Cov{C(f ). C(f ± )} = E{C(f ) C(f.)} - E{C(f.)} E{C(f.)} 


2 2 
2L a 


f+Kl-'HI'HI] 


(A- 


This equation can be evaluated with the help of Appendix 2, bearing 

2 

in mind that in this appendix the variance of A^ and B^ is yy instead 

- 2 
of a . 


r 


-14 2 

1.406 x 10 a A 


i=3 


Cov{C(f ), C(f.)} = i 


-242 i i 

6.2S0 x 10 o A , |i-j | = 1 


• (A- 


-1 4 2 

3.907 x 10 o A 


0 


i-j = 2 


> I i-j1 ^3 


As before, since non-overlapped data segments are independent, 
when 2 non-overlapped Hanned spectrum estimates are averaged, 


, . 4.2 

^ \ .406 x 10 1 2-A_ 


Cov{C(f ), C(f.)} 


6.250 x 10 


, 4 2 

-2 a A 


= < 


3.907 x 10° 


i=J 


i-j | = 1 


i-j] = 2 


» I i -j I ^ 3 


(A- 


3.9) 


3.10) 


3.11) 
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APPENDIX 4 

STABILITY OF ESTIMATED LOW RESOLUTION SPECTRA DERIVED FROM 
FIFTY PERCENT OVERLAPPED HIGHER RESOLUTION FFTs 
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i. 






In this appendix we will generalize the results from Appendix 2 
to the case where the high resolution spectrum estimate is derived from 
overlapped transforms. In particular, we will assume that the data has 
been divided into P segments where the last half of each segment over¬ 
laps the first half of the succeeding one. These segments are then 
windowed and transformed, and the resulting P spectra are averaged to 
form the high resolution spectrum estimate. The low resolution estimate 
is then derived from the high resolution estimate by means of the frequency 
averaging techniques. 

The following notation will be used in this appendix. Let x(t) be 
a realization of a stationary random process with power spectral density 
SCf). Let Cp(fj) be the discrete spectrum estimate based on the p tn 
segment of data. (Cp(f^) is formed as described in Chapter II by sampling 
the windowed data segment and doing a discrete transform with an FFT 
algorithm.) Then the high resolution estimate of S(fj) is 

P 

C(f.) = i-'V' C (f.) . (A-4.1 

r p p r 

p=i 

The frequency averaging techniques are employed to produce a low resolution 
estimate of S(f^) defined by 
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(A-4.2) 


e <v ■ £ a j C(f k-j ) 

j = -K 


where the are called the window coefficients for the frequency averaging 
on the high resolution spectrum. In order to compute the number of equiva¬ 
lent degrees of freedom v of C(f^), we will compute its mean and variance 
and make use of Eq. (2.26): 


(A-4.3) 


We first compute the mean: 

E{c( f k )}- E 

p=l j=-K 

To compute E[Cp(f^ j)] we must introduce further notation. 
w(t) be the data window, centered about the origin, and let 
length of each FFT in seconds. Then C p C f j) is the result of trans¬ 
forming the function x(t)w(t~^-p). To simplify matters, we will 
use continuous transform notation, the underlying assumption being 
that the results will not differ significantly if the sampling frequency 
is high enough to prevent corruption of the discrete transforms by 
aliasing. Then, from the autocorrelation theorem we have 


(A-4.4) 

Let 

L denote the 





2E 2 |£(f k ) 

Varic(f k ) 
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dudv, (A-4.S) 


E|c p (f k _j) I = ff e|x(u)x*(v)| w{u- w*(v-^-) e 


i2iT(u-v)f 


k-j 


where all integrals, unless otherwise indicated, are from -*> to °°, the 
asterisk denotes complex conjugate, and i=/^T. Also from the autocorrela¬ 
tion theorem. 


E{x(u)x*(v)} = J' S(x) e^ 27T ^ U V ^ T dx 


(A-4.6) 


Substituting into Eq. (A-4.5) ; 


= f S(T)|W(f k _.-x)| 2 dt 


f S(f R _j-x)|W(x)| 2 dx 


(A-4.7) 


where W(x) is the Fourier transform of the data window w(t). (W(x) is 

also referred to as the spectral window.) If we assume that the width 
of W(x) is such that S is approximately constant in a neighborhood of 
f^ j of that width, then we can make the approximation 


E { c P (f k-j ) }- S(f k-j ) / |ww|2dT 


(A-4.8) 



Further assuming that S is approximately constant over frequencies from 


f^ K to f k+K > we have from Eq. (A-4.4), 


={c(f k )} « SCf k )[/ |W(x )| 2 dx] 


_I\ 

LJ- K J 


(A-4.9) 
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Turning now to the computation of the variance of C(f^), we can 
expand in familiar fashion to get 


v ar{c ( f k )} = E{(c(f k )) 2 } - E 2 {c(f k )} 

Letting 


(A-4.10) 


C (f. ) = V" a - C (f. .) 

P ^ j p l k-J 

j = -K 


we can expand Eq. (A-4.10) further to yield 


r 

Hvv} 


^l[ E lwvvl- E {yv} E K (f i4] 


(A-4.11) 


Note that if the random process x is stationary, then 


Hyvl ■ Hv f k >} 


for all p and q. Similarly, 


E k ( v}- E { c , (f k>> 


for all p and q. Also note that for all p and q such that the corresponding 
segments of the windowed time function do not overlap 


={y f k » yv} ■ f-lyv} F {y<vi 










since C^(f k ) and C^(f k ) are independent for this case. Combining all of 
the above to simplify Eq. (A-4.11) we have 

Var{c ( f k )} ■Iv.rjytj)) * ( W f k>} ' E { C 'l< f k>} E { c_ 2‘ f k>l • ‘ A 

(For fractional overlaps less than or equal to 50%, there are 2(P-1) pairs 
of segments that overlap in Eq. (A-4.11) 

In order to evaluate Eq. (A-4.12) we will first compute the variance 
of Cj Cf k ): 



EE 

i=-K JL=-K 


V* 


E { c i (f k-j )c i (f k-t>> 


‘ K 

E a i E ( c i (f k-j 


j-K 


1 


(A-4.13) 


We have already computed E^C^f^)^ in Eq. (A-4.8). Using essentially 
similar techniques, we can show 


E { c i (f k-j> c i (f k-«>l' // s w s w 

X [|W(f k _. - Vi) | 2 |W(f k _ £ -x) | 2 

+ wc^.j-u) w*(f k j+ T) W(f k _^ + u) w*(f k . r t) 


4.12) 


+ W(f k _j-u) W*(f k _j-T) W(f k _ r x) W*(f k _ r p) dud! . (A - 4.14) 







Now assuming that the width of W is small in comparison with j and 
f k , we can neglect the second term in Eq. (A-4.14) since the windows 
will not overlap. 


E { C l( f k-j )C l( f k-W ff S(U)S(T)|W(f kj -u)| 2 |W(f k£ -T) | 2 dudT 

+ \f S(u)W(f k _ .-y)W* (f kJZ -y) djJ 


(A- 


Let n=f, .-p and o. = f .-f, Then, again using techniques already 

K-J J,K.K-JK-iC. 

demonstrated, 

f S(p) "(f k _ r ») ~-f S f f k .j~ n ^ w ( r 0 W* (n+o. z ) dn 

“ - S(f fc )/ W(n) V*(n+o. ) dn (A- 

Substituting back into Eq. (A-4.15), we can obtain 


; { C l (f k-j) ■ 5 2 <f k )[(/l»(T)| 2 d,)‘ 


*1/ 


W(t) W*(t+o. ) dx | 

3 »^ 


2 ] 


(A- 


Also necessary in the computation of Var C(fj.)is 


K K 

E { C ~l (f k> V f k>} " 2 2 V« E { C l (f k-j’ C 2 (f k-d} 

i=-K t=-K 


4.15) 


4.16) 


4.17) 


(A-4.18) 



where, letting k-j=n, k-£=m, 


E { C l (f n } C 2 (f m 3 } = //// E { x(u} X * (V) x(r) x * (s) } 


x w(u- E ) w*(v- w(r-L) w*(s-L) 


-i2n[f (u-v) + f (r-s)] 
x e dudvdrds 


(A-4.19) 


From [Ref. 6], we have 


e|x(u) x*(v) x(r) x* (s) \-ff S(y) S(x) 


x < e 


,i2n [y (u-v) + t (r-s) ] + e i27r[y(u-r) + t(v-s)] 


+ i2rr [y (u-s) + t (r-v) ] 


dydt 


(A-4.20) 


Substituting into Eq. (A-4.19) and integrating with respect to u, v, r, 
and s yields 

E { C ,< f „) ’ffsM S(,) j !«(f n -u) | 2 |W(f nl -T) | 2 


+ W(f -y) W*(f +t) W(f +y) W*(f - T ) e 
n n ' K m m 


— in(y+x) L 


+ W(f -y) W*(f -t) W(f -t) W*(f -y) e 
n n m v m 


- lit (y-i) L 


dydx 


(A-4.21) 
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5 


Again, the middle term of Eq. (A-4.21) can be neglected wherever both f 
and f are greater than the bandwidth of W. Thus, 




| f W*(t) W(r+f_-fJ e" 11,L dt 1 2 


m n' 


! ] 


(A-4.22) 


where we have again invoked the assumption that S does not vary 
significantly over the frequency range from f^^ to f^ + £. 

Combining Eqs. (A-4.8), (A-4.12), (A-4.13), (A-17), (A-4.18), and 
(A-4.22) yields 


v ar{c(f k )} • 


S2 (f k ) 


K K l 

E E a i a * «{iy»w i 

4 — is n _ v ' J 


dx 


j=-K a—K 


(A-4.23) 

The number of equivalent degrees of freedom of the frequency averaged 
spectrum can now be found from Eq. (A-4.3) with the aid of Eqs. (A-4.9) and 
(A-4.23). 


2 (P-1) 
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W(t) 


W*(t+o . ) 

J y * 


ilTtL 


dx 
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APPENDIX 5 

RESOLUTION REDUCTION USING COEFFICIENTS WITH 
UNMATCHED EQUIVALENT DEGREES OF FREEDOM 







In this appendix we will present examples illustrating the 
most commonly employed technique for resolution reduction. The procedure 
is to use frequency averaging with uniform coefficients, using the same 
number of coefficients as the factor by which resolution is to be 
reduced. For example, to reduce resolution by a factor of eight, eight 
uniform coefficients would be used. The difficulty with using this 
procedure is that neither the spectral window nor the statistics of the 
frequency averaged spectrum match those of a conventionally obtained 
spectrum. This can lead to a certain amount of misinterpretation of 
the results. 

In Figs. 29-34 we show, in the same format used in Chapter IV, 
examples of spectrum estimates resulting from a 4:1 reduction 
from a high resolution spectrum. The HSQ coefficients are the same ones 
used in the previous examples, while four uniform coefficients were used 
to produce the spectrum shown in the bottom plot of these figures. From 
the appearance of these plots, we could conjecture that the spectral 
window associated with the four uniform coefficients is somewhat narrower 
than the Hanning window used to produce the conventional low resolution 
estimate. The spectral window comparisons for this type of resolution 
reduction are shown in Figs. 55-37 for 2:1, 4:1, and 8:1 reduction ratios. 
An examination of Fig. 36 shows that the window for this 4:1 example is 
indeed narrower than the Hanning window. Hence, our resolution lias not 
actually decreased by a factor of four, as desired. 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


FIGURE 29 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, LxF = 32.000 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 
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FIGURE 30 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, L x F = 32.100 
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HIGH RESOLUTION SPECTRUM 



APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


FIGURE 31 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, LxF -32.200 
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APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


FIGURE 32 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, L x F - 32.300 
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APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


FIGURE 33 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, L x F - 32.400 
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APPROXIMATION OF LR SPECTRUM WITH HSQ COEFFICIENT 



APPROXIMATION OF LR SPECTRUM WITH UNIFORM COEFFICIENT 


FIGURE 34 

COMPARISON OF POWER SPECTRA 
4:1 REDUCTION, L x F - 32.600 aEwmk 
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FIGURE 35 

SPECTRAL WINDOW COMPARISONS, 2:1 REDUCTION 
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FIGURE 36 

SPECTRAL WINDOW COMPARISONS, 4:1 REDUCTION 














FIGURE 37 

SPECTRAL WINDOW COMPARISONS, 8:1 REDUCTION 
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Examination of the spectra shown in Figs. 29-34 also suggests 


that the uniform coefficient frequency averaged estimate lacks the 
stability of the other low resolution estimates. Computation shows 
that the conventional low resolution estimate has 89 equivalent degrees 
of freedom, while the four uniform coefficients produce only 48. This 
corresponds to an increase in minimum detectable signal level of approxi¬ 
mately 1.5 dB. However, over the frequency range of interest, the 
uniform coefficient window averages approximately 1.2 dB greater than 
the Hanning window. Hence the net degradation in minimum detectable 
signal level due to the altered spectral window, neglecting equivalent 
noise bandwidth changes, is approximately 0.3 dB. 
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